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A detailed knowledge of the gap structure for the Fe-pnictide superconductors is still rather rudimentary, with 
several conflicting reports of either nodes, deep gap minima, or fully isotropic gaps on the Fermi surface sheets, 
both in the k x — k y plane and along the c— axis. In this paper we present considerations for electronic Raman 
scattering which can help clarify the gap structure and topology using different light scattering geometries. 
Using density functional calculations for the Raman vertices, it is shown that the location of the gap minima 
may occur on loops stretching over a portion of the c-axis in Ba(Feo.94Coo.o6)2As2- 

PACS numbers: 74.25.nd,74.70.Xa,74.20.Pq,71.15.Mb 



Since the discovery of high temperature superconductivity 
in the iron pnictides, identifying the structure of the super- 
conducting order parameter has remained a tantalizing and 
elusive problem. Knowledge of the momentum dependence 
is a major step towards identifying the interactions that drive 
the superconductivity. A variety of experimental techniques 
probing the low energy properties of the pnictides in the su- 
perconducting state have given contradictory information on 
the angular dependence of the superconducting gap around 
the Fermi surfaces (FS). Angle-resolved photoemission on 
NdFeAsOo.gFo.i 1 , K-doped 2 and Co-doped 3 BaFe2As2 have 
identified rather isotropic superconducting gaps on both the 
hole and electron pockets in the Brillouin zone (BZ). On 
the other hand, transport and thermodynamic measurements, 
which do not have angular resolution, have indicated the pres- 
ence of nodes in some materials and small but finite gaps 
in others 4 . Recently, angle-resolved thermal conductivity in 
BaFe2(As,P)2 have indicated large gap variations around the 
Fermi surface, with gap nodes or deep minima located at 45° 
with respect to the principal BZ axes 5 . 

Electronic Raman scattering experiments indicated very 
small gaps located on the electron Fermi surfaces in 
Ba(Feo.94Coo.o6)2As2 6 . However, due to the relative insen- 
sitivity of the Raman vertices to momentum in the canoni- 
cal tight binding (TB) scheme the location of the nodes could 
not be determined uniquely. The situations is additionally ob- 
scured by the fact that the electronic structure of the 122 mate- 
rials is inherently more 3D than that of the 1111 pnictides, and 
simplified 2D tight binding models can be rather misleading 7 . 
Moreover, the crystallographic symmetry of the 122 structure 
is body-centered tetragonal (bet), IA/mmm, therefore folding 
down the Fermi surfaces from the one-Fe Brillouin zone to the 
two-Fe zone is less trivial than in the 1111 structure. Namely, 
while in the 1111 structure the electronic pockets at the X 
point and Y point in the unfolded zone are related by a simple 
Tt/2 rotation, so that after folding they form e-pockets with the 
full 4-fold symmetry, in the 122 structure the symmetry oper- 
ation connecting the two unfolded pockets is the same rotation 
plus a shift along k z by n/c. Thus for all k z but k z = 7t/2c the 
folded FS cuts have only a 2-fold symmetry. 



Incidentally, for the same reason the standard method of 
constructing a TB Hamiltonian by transforming the Bloch 
functions into the Wannier functions leads to a TB model that 
attempts to describe the outer barrel as a single band in the 
unfolded BZ 7 , while in reality it is formed by anticrossing the 
X and the Y pockets. In other words, such TB models feature 
unphysically and artificially large interlayer hoppings (some 
interiayer hoppings in Ref. 7 appear to be more than 50% 
of the nearest-neighbor in-plane hopping), while in reality the 
scale of the k z dispersion is set by the hybridization gaps, and 
by the difference in size (but not orientation) between the up- 
per and the lower panel in the Fig. 4 of Ref. 7 . 

For this reason, in the 122 systems extracting quantitative 
information from the experiment using an infolded 2D tight 
binding model becomes questionable, and working in a folded 
3D TB model is not much simpler than straight off first prin- 
ciple calculations. Therefore, in order to gain more accurate 
information about the location of gap minima on the Fermi 
surface we have performed all-electron band structure calcu- 
lations, as described in Ref. 8 and numerically computed the 
relevant Raman vertices. We have confirmed the result ob- 
tained qualitatively by Muschler et al. 6 that two of the three 
principal Raman symmetries, A\ g and Big (all notations in the 
crystallographic unit cell), probe, respectively, an average gap 
and the gap on the electron Fermi surface, and that for the third 
symmetry, B\ g , the matrix elements are mostly small on all the 
FSs. However, we found that on the outer electronic FS there 
exist "hot spots" where this matrix element is sizeable, a fact 
that could not be seen in unfolded TB calculations. Therefore 
we can not only distinguish between the hole and the electron 
FSs in terms of the gap minima, but can also place them at or 
near these hot spots (which appear near the [1 10] direction). 

First, we discuss the so-called "effective mass approxima- 
tion" as applied to this compound. The inelastic light (Raman) 
scattering cross section in the limit of small momentum trans- 
fers is given in terms of correlation functions of an effective 
charge density 9 , p = Y,n.k,° T«( k ' ®i, ®s)c\ k a c„ ik!CT , where 
c n k a ( c n,k,a) creates (removes) an electron with momentum k 
and spin a from band n. In principle, y depends both on light 
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FIG. 1 : (Color online) Raman response for xx and xy polarization for 
three different excitation lines as indicated. All spectra have the Bose 
factor removed. Apart from correcting for the instrumental response, 
the overall intensity is as measured and has not been shifted. 



scattering polarization geometries e'' s ' and the frequencies of 
incoming (scattered) light (o'( s K In the limit where (D'W are 
much smaller than any relevant interband transition frequency, 
the scattering amplitude simplifies into the well-known effec- 
tive mass tensor 



y„(k,(» J >^0) 



a,P=x,y,z 



' a dk n dk/p- 



(1) 



For tetragonal materials, three different scattering geometries 
can be classified according to the point group transformation 
properties. These are given as 

A lg : y„(k) = [d 2 e„(k)/dk 2 x + d 2 e n (k)/dk 2 }/2 
B lg : y„(k) = [d 2 e n (k)/dkj - d 2 e n (k)/dk 2 }/2 
B 2g : Y„(k) = d 2 e n (k)/dk x dk y . (2) 

Thus the curvature of the bands and the light scattering po- 
larization orientations determine which carriers are involved 
in light scattering in different bands and regions of the BZ. 
While not having the definitive angular resolution as photoe- 
mission to probe single particle properties, Raman scattering 
can selectively project effective charge density fluctuations in 
different regions of the BZ, uniquely providing k— projected 
dynamical two-particle properties. 

In most systems, such as the cuprates and pnictides, inci- 
dent light energies are typically high enough to directly ex- 
cite interband charge transfer and thus the effective mass ap- 
proximation may be not valid for these systems. When this 
happens, strong resonance effects occur and the intensity of 
the electronic Raman scattering varies greatly with the inci- 
dent photon energy, as for instance in cuprates. 10 However, 
as shown in Fig. 1 in the pnictides the resonance effects are 
weak and thus off-resonance conditions for the scattering am- 
plitude applies. Therefore we assume that the effective mass 
approach holds for the pnictides and examine the full band 
and momentum dependence of the scattering amplitudes as a 
way of pinpointing the gap extrema structure. 



In the superconducting state, the Raman response for B\ g 
and Bi g symmetries is given as a weighted average of the 
Tsuneto-Maki function A„ over the FS for each band n , 
%y s y{Q.) = £„. k y 2 (k)A„(k,£2), with the Raman vertices y 
given in Eq. 2. The Tsuneto function is 

E n (k) 4|A„(k) \ 2 /E n (k) 



A„(k,£2) = tanh 



2k B T 4£ 2 (k)-(M2 + ;<5) 2 



(3) 



where A„(k) is the superconducting gap in band n, having 
quasiparticle energy dispersion £„(k), and E 2 ik) = £ 2 (k) + 
A 2 (k) . Enforce particle number conservation, one gets 

X(0) = Xy.y(&) - Xy.\ («)Zi, y («)/Zi,i («)> (4) 

with^i (£2) =Xi,y,{&) =L,, k y,(k)A„(k,£2) and^ u (£2) = 
k A,,(k,£2). Therefore, while for both B\ g and B 2g each 
contribution to the Raman vertex can be separated into con- 
tributions from single bands, the A\ g contributions mix mass 
fluctuations in each band and between bands. 

In Ref. 6, the interplay of the energy gap and the symmetry 
dependence of the Raman vertices was considered, where it 
was argued that B2 g polarizations largely sampled the energy 
gap on the electron bands, while A \ g was more sensitive to the 
hole bands in the pnictides. Here we explore the form for the 
Raman vertices in more detail using LDA. 

The Raman vertices shown in Eq. (2) are directly com- 
putable from the one-electron band structure. We have per- 
formed full-potential LAPW calculations using the WIEN2k 
package and an experimental crystal structure, including the 
spin-orbit interaction. The doping is accounted for in the 
virtual crystal approximation. The derivatives of the eigen- 
states were calculated numerically using a 31x31x31 fc-point 
mesh in the bet BZ. We neglect any momentum interplay be- 
tween the Raman vertices and the superconducting gaps in 
the Tsuneto-Maki function, and examine just the character 
of the Raman vertices on each band for each polarization 
geometry. 12 The calculated quantities are derived from eval- 
uating Eq. 4 as 



x; c ; r =j:^(m^)]- 



{L n . k Uk)8[e n (k)]y 
E„,kS[£„(k)] 



(5) 



This can be simplified by writing y„(k) = Ay„(k) + %, with 
?n = Lk7n(k)<5[£„(k)]/£ k <5[£„(k)] the average of the Raman 
vertex for band n over Fermi surface sheet n. The response 
can be separated into an intraband piece, which is the only 
contribution for B\ g and Bi g , 



Xintraband = £[Ay„(k)] 2 (5[e„(k)], 
n.k 



(6) 



and an interband contribution only appearing for A\ g channels 
UnNn? 2 — (LnN„7n) 2 /LnMz- The last part can be written as a 
band-dependent sum: 



Xinterband ,t Yn / ,^m(Yn Ym) — T^. Xn, interband t 

n t m n 



(7) 



with N t = £„iV„ the total DOS. The results are shown in Table 
I and are defined as follows: 

JV n = £S[e n (k)] 

k 

7n,ap = £[d 2 5,(k)/dfc a dty]5[E„(k)]. 

k 

£ [d 2 e„ (k)/dk a dkp][d 2 e n (k)/dk Y dk 5 \ 8 [e n (k)] , 

k 

such that the B\ g polarizations select (y„, xxxx — Yn, XX yy) /2, Big 
selects y n . X y X y The unscreened intraband part for A\ g is given 
by {Ynjcxxx + Yn xtyy)/2> which is reduced with screening to 

(y«^xa 4" Yn,xxyy) /2 N n f ■ 
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TABLE I: Values of the Raman vertices and related quantities in- 
cluding spin-orbit coupling around the hole Qil,h2,h3) and electron 
{el, el) Fermi surfaces, on a per cell basis, in Ry units (energy in 
Ry, distance in Bohr radii, momentum in inverse Bohr radii). Here 
hi ...Kije\...e2 denote the outermost.. .innermost hole/electron barrel 
Fermi surface, centered at TIM in the Brillouin zone, as indicated in 
Fig. 2. 

The LDA results shown in Table I indicate that both B\ g 
and B2 g polarizations give vertices which are dominated by 
contributions from the electron FS, in agreement with the as- 
sessment of Ref. 6. Moreover we note that B\ g vertices are 
specifically large only on a single electron FS sheet, given 
by the outer FS pocket centered at M. For the A\ g contri- 
bution, while the unscreened Raman vertex y xxxx + y xxyy has 
large contributions from the hole band hi and electron band 
el, backflow largely cancels the intraband contribution, leav- 
ing a substantial interband contribution deriving largely from 
these same bands. These results support the conjecture based 
on symmetry in Ref. 6 that the hole bands would be probed 
for A\ g configurations more so than for B\ g and Big- 

The LDA results indicate that the A\ g interband portion in- 
volves also a strong contribution from scattering between the 
electron band el and the hole band hi, coming from different 
signs of the Raman vertices in the two bands. As discussed 
in Ref. 1 1 the divergence at twice the gap edge which may 
appear in crossed polarization orientations (such as B\ g 2 g ) is 
removed in single band superconductors for Ai g symmetries 
due to the gauge invariant backflow which screens uniform 
charge displacements over many unit cells. In multiband sys- 
tems, such as the pnictides, the superconducting response for 
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FIG. 2: (Color online) Left: Raman vertices for B\ g symmetry, eval- 
uated around each Fermi surface sheet projected onto the x — y plane, 
that is, looking down the ^-direction. The parts of the Fermi sur- 
face where the value of B\ g Raman vertex is larger than the average 
are shaded red. The gold ellipses mark the area where the deep gap 
minima or nodes are expected. Right top: The same, viewed along 
the M — r line, as shown by the arrow in the left panel; only the 
outer electron FS pockets are shown. Right bottom: The same view, 
but now the parts where the Fermi velocity is larger than 60% of its 
maximal value are red. 

A\ g symmetries may be written as 11 

x s a% (A) = E (Yn (k) - ym 2 K (k, a) , (8) 

n,k 

with X„ given by Eq. (3). The Tsuneto-Maki function A„ pos- 
sesses divergences when Q. = 2 |A„(k)|. y(il) is another way 
of representing the frequency-dependent backflow in the su- 
perconducting state for A\ g symmetries: 

y(O) =£y„(k)^(k,n)/£\(k ! n). (9) 

«,k /7,k 

Each separate divergence occuring at the gap edge in each 
band in A„(k,£2) weighs exactly those momentum points in 
y„ (k) for which the divergence occurs, and via Eq. 8 the di- 
vergence is cancelled. This holds for each divergence coming 
from each band n, and therefore we conclude that no diver- 
gence may appear in Ai g symmetries for multiband systems. 
This conclusion holds only for Cooper pair creation via light 
scattering but does not consider additional divergences com- 
ing from charge or spin collective modes. 13 ' 14 

This predicted absence of an A \ g divergence is in general 
agreement in the observed data on the pnictides, where no 
clear 2A scattering peak is seen for A\ g scattering geome- 
tries. While sharp collective modes may appear due to the 
oscillation of pairing amplitudes in the different bands for A\ g 
symmetries 15 , no clear indication of such a mode has been 
found. Similarly, although collective modes have been pre- 
dicted to also occur in Big symmetries due to final state exci- 
tons having orthogonal symmetry to the superconducting pair 
state 16 , it is still an open question whether these modes can be 
uncovered in the data. 

The fact that in the experiment hardly any effect of super- 
conductivity is observed in the B\ g polarization, while the cal- 
culations predict a noticeable coupling with the outer elec- 
tronic FS, suggests that this coupling is not uniform, but is 
small wherever the gap is sizeable, and where the coupling 
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is large the gap is small. In order to test this hypothesis, we 
have plotted the effective mass asymmetry, appearing in B\ g 
polarizations as \d 2 e n (k) /dk^ — d 2 e n (k)/dk 2 \, in Fig. 2. We 
observe that, first, this difference is small on all FSs except the 
outer e-barrel, and, second, it is also small on the most of that 
FS as well, except on a small spot located at an angle varying 
from 10 to 15", depending on k z . As Fig. 2 shows, this is con- 
sistent with gap nodes or deep gap minima, forming elliptical 
loops on the outer electron FS centered around the (£, , £ , 0) 
points, where the F — M line intersects the electron pockets 
(indicated by the ellipse in Fig. 2). This can be compared with 
model calculations by Graser et al. 7 who found that the pair- 
ing interaction is strong in the Z plane (k z = n/c) plane and 
weaker in the F plane. In their calculations this led to a nodal 
loops similar to ours, but located at the hole FS. One should 
be, however, careful, when making a quantitative, as opposed 
to qualitative, comparison, with these calculations, because, 
as discussed above, Graser et al. worked in an unfolded BZ 
wher the folding-induced 3D dispersion was emulated through 
enhance k z hopping parameters. 

It is also consistent with the thermal conductivity data 5 
(collected on another 122 material), which indicate that there 
are gap nodes present in the system, yet the specific heat does 
not seem to be affected by them. These data have been in- 
terpreted in terms of ultraheavy hole bands that contribute at 
least 80% to the DOS, and ultralight electron bands that dom- 
inate the transport. However, the former assumption seems 
a bit too strong to be realistic. Yet if in BaFe 2 (As,P) 2 the 
gap structure resembles the one that we propose here for 
BaFei.8gCoo.i2As2, so that the gap nodes exist only for a finite 
range of k z , the rest of the electronic bands, fully gapped, add 
their contribution to the hole band and together provide 80% 
or more of the total specific heat. In fact, incidentally, the 
Fermi velocity is the largest close to the "hot spots" that we 
have identified as locations of gap nodes, in agreement with 
the idea that the transport is dominated by near-nodal parts of 
the FS (Fig. 2). 

One can ask a question, why most of the Fermi surface, 
except for relatively small "hot spots", shows such weak cou- 
pling in the B\ g channel. This is actually easy to address. Let 
us start with a simple model, a 2D parabolic electronic band. 
By symmetry, the principal axes of the corresponding ellipti- 
cal FS are along (11) and (11). For the moment, we will use 
the center of this FS as the origin and the principal axes as 
coordinates. Then £12 = Mi,2&* + M2,i&y- The Raman vertex 
in question appears like a Big vertex in this basis and is pro- 
portional to d 2 e„(k x ,k y ) I dk x dk y = 0. Thus in this model this 
band does not couple in the B\ g channel at all. 

However, this is not necessarily true if one takes into ac- 
count the downfolding of the electron pockets in the crystal- 
lographic unit cell. After the downfolding the two FSs simply 
intersect. In our model, as well as in the real 1 1 and 1111 
crystallography, they cross along the k x = and k y = di- 



rections, and, because they have different parities, there is no 
hybridization between them and the newly formed degener- 
acy is not lifted (although it would be lifted by the spin-orbit 
coupling). On the other hand, in the 122 symmetry the inter- 
sections occur at general points in the Brillouin zone so that no 
symmetry prevents them from hybridization, thus forming an 
inner barrel and an outer barrel (el and el) even without the 
spin-orbit. To gain more insight into this problem we intro- 
duce hybridization into our model, and write the Hamiltonian 
near the anticrossing point as H\\ = E\ (k),H22 = £2(6) ,#12 = 
H21 = V, where V is a hybridization matrix element. The in- 
verse masses of the resulting two bands are 



d 2 e ± 



dk x dk y 



SV 2 \(m-H2) 2 k X ky 



[( J U 1 - Al2 ) 2 fe 2 -L2)2+4y2]3/2- 



(10) 



This is essentially an approximate expression of a delta func- 
tion centered at k x ± k y = with the width V and strength 
\H\ — jU 2 |. Thus in this simple model we have vanishing cou- 
pling everywhere except hot spots where the bands would 
cross if the As sublattice would not have a lower symmetry 
than the Fe one. 

This simplified picture only qualitatively applies to Bal22. 
The bands there are not exactly parabolic, and, most impor- 
tantly they cross not at the k x , k y = lines, but at general points 
that can be considerably removed from these line, depending 
on k z . Yet our model demonstrates qualitatively how the "hot 
spots" for this polarization appear on the FS, and why they are 
so small in size. 

In summary, in this paper we have paired experimental data 
on the electronic Raman scattering in the superconducting 
state of Ba(Feo.94Coo.o6)2As 2 with the first principle calcu- 
lations of the effective mass fluctuations. We found that a 2D 
TB model in the unfolded Brillouin zone is not sufficient to 
capture all the mass fluctuation. The reason is that because 
of the body-centered symmetry of the corresponding crystal 
structure downfolding the electron bands results in hybridiza- 
tions and substantial band anticrossings, which in turns pro- 
duces large mass fluctuation in a small part of the Fermi sur- 
face. Our results suggest that closed nodal loops (rather than 
vertical or horizontal nodal lines) exist in this compound near 
the intersection of the outer electron Fermi surface with the 
T— M line. It is likely that the size and the very existence of 
these loops is doping-dependent, but this is beyond the scope 
of this paper. 
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